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We examine the constraints on final state radiation from Weakly Interacting Massive Particle 
(WIMP) dark matter candidates annihilating into various standard model final states, as imposed 
by the measurement of the isotropic diffuse gamma-ray background by the Large Area Telescope 
aboard the Fermi Gamma-Ray Space Telescope. The expected isotropic diffuse signal from dark 
matter annihilation has contributions from the local Milky Way (MW) as well as from extragalactic 
dark matter. The signal from the MW is very insensitive to the adopted dark matter profile of the 
halos, and dominates the signal from extragalactic halos, which is sensitive to the low mass cut-off 
of the halo mass function. We adopt a conservative model for both the low halo mass survival 
cut-off and the substructure boost factor of the Galactic and extragalactic components, and only 
consider the primary final state radiation. This provides robust constraints which reach the thermal 
production cross-section for low mass WIMPs annihilating into hadronic modes. We also reanalyze 
limits from HESS observations of the Galactic Ridge region using a conservative model for the 
dark matter halo profile. When combined with the HESS constraint, the isotropic diffuse spectrum 
rules out all interpretations of the PAMELA positron excess based on dark matter annihilation into 
two lepton final states. Annihilation into four leptons through new intermediate states, although 
constrained by the data, is not excluded. 



I. INTRODUCTION 

The identity of the dark matter has remained a fun- 
damental unsolved problem in cosmology and particle 
physics for nearly eighty years PQ. The likelihood of 
new physics near the electroweak scale suggests a par- 
ticle dark matter candidate of mass ^100 GeV. A weak- 
scale interaction strength for such a particle can nat- 
urally produce a relic abundance at the observed dark 
matter density. Such a weakly-interacting massive par- 
ticle (WIMP) arises in supersymmetric extensions of the 
standard model as well as other beyond the standard 
model extensions. For a review of particle dark matter 
candidates, see, e.g., Refs. [2 [3]. 

Thermal production of dark matter prefers a scale of 
the dark matter annihilation cross-section at {u^v) w 
3 x 1CP 26 cm 3 s -1 . This scale predicts a potentially de- 
tectable gamma-ray photon signal from annihilations in 
local and cosmological dark matter structures, e.g., see 
Ref. [4]. Annihilation channels directly to photon pairs 
produce a distinctive 7-ray spectral line, while channels 
to charged standard model particles produce an associ- 
ated continuum emission of 7-rays from bremsstrahlung 
radiation and in hadronization via tt° — > 77. The Large 
Area Telescope (LAT) aboard the recently launched 
Fermi Gamma-Ray Space Telescope has significant sen- 
sitivity to such 7-ray radiation from dark matter annihi- 
lation. The sensitivity of Fermi-LAT to a an annihilation 
signal from the Galactic center (GC), Galactic satellites 



and the isotropic diffuse signal was reviewed in Baltz et 
al. [5]. 

Recent Fermi-LAT observations have led to constraints 
on dark matter annihilation, including observations of 
dwarf galaxies analyzed by the Fermi-LAT Collabora- 
tion [5J [7] , as well as analyses of public data of the total 
sky flux [H], and public data from selected portions of 
the sky [9]. Most recently, the Fermi-LAT collaboration 
released the spectrum from the isotropic diffuse back- 
ground in 7-rays [TO]. Here, we use the Fermi-LAT Col- 
laboration derived diffuse-photon spectrum in this work 
to place constraints on dark matter annihilation chan- 
nels. Depending on the background estimation, the rela- 
tive sensitivity of these methods varies. 

In this work, we constrain the partial cross section for 
dark matter annihilation to various particle anti-particle 
standard model final states using the diffuse isotropic 
7-ray photon spectrum from the Fermi-LAT 10J . The 
Fermi-LAT isotropic diffuse spectrum uses enhanced 
event simulation modeling that improves the background 
rejection in the observations by a factor of 1.3-10. There- 
fore, we expect the spectrum to be a stringent constraint 
on dark matter annihilation. The errors on the estimate 
of the diffuse isotropic background come from system- 
atic errors in modeling foreground Galactic diffuse 7-ray 
emission. The isotropic diffuse background could include 
a dark matter annihilation signal, yet also has contri- 
butions from unresolved active galactic nuclei, starburst 
galaxies and 7-ray bursts, as well as any truly diffuse 
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components in large scale structure. In total, the model- 
ing of the Galactic diffuse emission, as well as the extra- 
galactic 7-ray emission processes, are the ultimate obser- 
vational limits to the sensitivity of the isotropic diffuse 
background to a dark matter signal. 

Since Cold Dark Matter (CDM) particle candidates 
such as the WIMPs considered here produce cusped ha- 
los [H[TT], in applying our constraints, we use conserva- 
tive models for the expected dark matter profile density 
distribution which are consistent with the cusped profiles 
produced by such CDM candidates. We also apply this 
requirement to update limits from existing constraints 
from observations of the GC ridge by HESS. 

The Fermi-LAT collaboration performed a similar 
analysis of the constraints from the diffuse background 
observations |12j . That work did not include the mini- 
mal Galactic contribution that is necessarily present in 
the observation, and is therefore less stringent than the 
limits presented here. Previous work along these lines 
prior to Fermi-LAT data include limits from older 7-ray 
observations with a subset of standard model channels: 
Bell & Jacques [13] examined constraints on the branch- 
ing ratio to e+/e~ with the EGRET, CGRO, HESS and 
CELESTE data, and Mack et al. [T3] applied constraints 
on lines in the case of the two-photon channel with a sim- 
ilar set of data. With the dawn of the Fermi-LAT data, 
constraints from all-sky and diffuse flux limits have been 
made with public data. Cirelli et al. jS] use the total flux 
from specific sky regions to constrain dark matter anni- 
hilations using conservative models for the dark matter 
halo profile. They include 7-rays arising from inverse- 
Compton scattering (ICS) of annihilation products on 
cosmic microwave background and inter-stellar photons. 
In a similar analysis, Pappucci & Strumia [8] find limits 
to dark matter annihilation modes from a full-sky anal- 
ysis of the public Fermi-LAT data, including final state 
radiation limits as well as ICS, with no reduction of the 
sky flux due to foreground galaxy and point source re- 
moval. Therefore, their constraints are not as stringent as 
ours or those in Ref . [5] . Ref. [5] also overstates MW halo 
profile dependencies by adopting isothermal halo profile 
models, which, as we discuss below, are inconsistent with 
the CDM structure formation paradigm for a WIMP can- 
didate. 

There has been considerable interest in the possibility 
of dark matter annihilation as being the source of the 
excess cosmic ray positron fraction at ^10-100 GeV ob- 
served by HEAT 1 [15], AMS-01 2 [TB] and PAMELA 3 Q7], 
where e + /e~ pairs are produced directly or indirectly 
in a dark matter particle pair annihilation cascade [181 
[T9] . In addition, features in the higher-energy 10 2 - 
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3 Payload for Antimatter Matter Exploration and Light-nuclei As- 
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10 3 GeV positron spectrum seen by ATIC 4 [2U] and 
Fermi-LAT [21] are also consistent with the dark mat- 
ter annihilation interpretations of in the lower energy 
positron excess data [221423] . In order to achieve the 
dark matter annihilation rate required for these e + /e _ 
signals that may be consistent with the expected ther- 
mal production cross section, and to avoid an excess in 
anti-proton observations, the annihilation rate can be en- 
hanced through a low-energy Sommerfeld enhancement, 
and limited to leptonic modes with a <1 GeV dark- force 
carrying particle [25- 28J . Such an enhanced cross-section 
from a new force is in tension with detailed calculations 
of the relic abundance of the dark matter, so that such 
a candidate may not contribute to all of the dark mat- 
ter [251-31 . and are also constrained by nonthermal dis- 
tortions of the CMB [29 and asphericity observed in dark 
matter halos [30] . 

These models are also constrained by corresponding 
Fermi-LAT 7-ray observations of dwarf galaxies [7], the 
total sky flux [5], portions of the sky [S], Galactic ra- 
dio synchrotron emission |32j , the neutrino flux from the 
GC as observed by SuperKAMIOKANDE [21], and at- 
mospheric Cerenkov observations of 7-rays from dwarf 
galaxies |55] . Hiitsi et al. [M] have used diffuse 7-ray ob- 
servations by EGRET to constrain dark matter models 
that can explain the PAMELA data, but as we explain 
below, with more optimistic assumptions for the low 
mass halo cutoff and extragalactic signal. The PAMELA, 
HEAT, and Fermi e + /e~ signals are also consistent with 
high-energy e + /e~~ emission from pulsars and supernova 
remnants 35 38 . Here, we show how the dark mat- 
ter annihilation interpretation of theses signals is in con- 
flict with the observed isotropic diffuse flux spectrum 
of Fermi-LAT in combination with other constraints for 
two-body standard model particle final states, and con- 
strains scalar or vector boson mediated four-lepton final 
states. 



II. ANNIHILATION SIGNAL 

A robust calculation of the expected final state radi- 
ation from dark matter annihilation requires accurate 
quantification of the dark matter source as well as the 
products in the final state 7-ray radiation chain. The 
observable flux of photons per solid angle for a specific 
annihilation channel with cross section (<jav) is 

_ (a A v) Jah 1 dN 7 
dE ~ 2 Jo AfUsTO 2 dE' U 

where J7ao /Jo is the normalized integral of mass density 
squared of the dark matter in the field of view, dN 1 /dE is 



4 Advanced Thin Ionization Calorimeter, 

http://atic.phys.lsu.edu/aticweb/ 
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FIG. 1: Shown is the observed isotropic diffuse background 
spectrum observed by Fermi-LAT (points with errors), and 
representative models of the annihilation spectrum of the 
Galactic (upper) and extragalactic (lower) contributions from 
the channel \X for m x = 10 GeV (solid), and for 

m x = 100 GeV (dashed). Both cases have the annihilation 
cross section (uad) = 3 x 10~ 26 cm 3 s _1 . 



the 7-ray spectrum per annihilation, m x is the dark mat- 
ter particle mass, and ASl b s is the observational solid an- 
gle in steradians. This relationship is dependent on the 
astrophysically-inferred dark matter density contribution 
in the field of view of an observation J7aq, and the parti- 
cle physics of the expected photon flux for a specific par- 
ticle candidate annihilation mode, dN 1 /dE. Note that 
dQ-JdE is the total photon number flux per unit energy 
per unit steradian for a full sky observation, and when 
compared to the total photon count of the Fermi-LAT 
observation with |6| > 10°, must be scaled to the field of 
view of that observation, A£! b s = 10.4 sr. 5 The Fermi- 
LAT diffuse isotropic background spectrum observation 
as well as examples of the isotropic diffuse signal for both 
the Galactic and extragalactic contributions are shown in 
Fig. [I] We describe our methods in detail below. 



A. Particle Annihilation Event Modeling 

For dark matter annihilations to two-body standard 
model final states, in order to obtain the average number 
of photons above a given energy per annihilation event 
and the differential number of photons at a given energy 
per annihilation event dN^/dE, we use Pythia 6.4 [39] to 
simulate both photon radiation off of charged particles as 
well as decays of particles such as the tt°. Specifically, we 
run Pythia to simulate an e + -e~ collision at a center of 
mass energy of 2m x through a Z' to a final state that cor- 
responds to the annihilation products of the dark matter. 



5 An earlier version of this paper incorrectly omitted this geomet- 
rical factor. We thank Julie McEnery, Michael Gustafsson, and 
the Fermi Collaboration for helping resolve this issue. 



For certain final states such as gg we use an s-channel h 
rather than a 2', as the corresponding decay process is 
already implemented in Pythia and does not need to be 
added by hand. For the hh final state we add a new decay 
channel for the Z'. We switch off initial state radiation 
such that all photons are emitted either radiatively off 
the final state particles or the decays of unstable parti- 
cles such as mesons. We use the default Pythia cutoff 
values for photon emission: quarks are assumed not to 
radiate below a GeV, and leptons below 100 keV. We 
turn on the decays of particles which are not decayed 
with the default Pythia settings, such as muons, charged 
pions and kaons. Using a large sample of events for each 
final state and each value of m x , the number of photons 
in the final state above a given energy are counted and 
averaged over the number of events, yielding the average 
number of photons above a given energy per annihilation 
event. 

For dark matter annihilations to final states such as 
4e and 4/j, through intermediate light scalars, we utilize 
the formulas given in Appendix A of Ref. [40]. For the 
four-e final state, all photons originate from final state 
radiation off of the electrons. For the four-/i final state, 
photons can originate from final state radiation off of the 
muons, as well as from final state radiation off of electrons 
produced by the decay of the muons. Since dark scalars 
are by assumption light, the energy range over which 
photons can be radiated is narrow in the rest frame of 
the scalar. While in the center of mass frame of the 
dark matter annihilation some of these photons can still 
carry a significant fraction of the available energy, the 
average number of hard photons per annihilation event is 
significantly reduced in comparison to direct annihilation 
to a two-body standard model final state. 

B. Cold Dark Matter Halo Models 

Thermally-produced dark matter particle candidates 
detectable in 7-rays are CDM particles whose primor- 
dial velocities are only determined by the initial gravi- 
tational perturbations. We therefore only consider dark 
matter halo models that arise in CDM cosmologies. In 
particular, density profiles of approximately the canoni- 
cal Navarro- White-Frenk (NFW) profile are expected for 
the Milky Way (MW) dark matter halo [TTJ , with profiles 
of the form 

pNPw(r) = p s ( ~\ + > ( 2 ) 

where p s is the characteristic density and r s is the scale 
radius. Halos in simulations exhibit scatter about this 
profile [UJ . We employ the NFW profile in order to allow 
for comparison to other work. However, higher resolution 
simulations of CDM halo formation have revealed a soft- 
ening of the profile power law with decreasing radius that 
deviates slightly but significantly from NFW, with a log- 
arithmic slope that decreases with radius [HI 02] ■ Called 



4 



Thermal Dark Matter 



100 1000 10* 

m y [GeV] 




HESS GR 
Thermal Dark Matter 







1 


XX -* uu ~- 

i i 


10 


100 1000 10 


- , 1 II III 1 


m x [GeV] 




HESS GR _ w i n- 24 _ 



1000 
[GeV] 




FIG. 2: Constraints on the partial cross section of annihilation of dark matter to quarks, leptons, gauge bosons and Higgs, for 
rrih — 120 GeV. The regions are labeled by their corresponding constraining observations as described in the text: "Fermi-LAT 
Diffuse" from the Fermi-LAT isotropic diffuse 7-ray background. The regions labeled "HESS GR" are for three different cases, 
solid cyan using the conservative Einasto profile as described in the text, hashed cyan from the most stringent case, a low 
mass high concentration NFW profile, and the canonical p Q = 0.3 GeV cm -3 , Rq = 8.5 kpc NFW profile as a dotted line for 
reference. All constraints are at 95% CL. 
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an Einasto profile, it is of the form 

2 ffr^ aE 

OLE \\T S 



PEinastoM = Ps exp 



- 1 



(3) 



where ctE is fit by simulations to be roughly ue ~ 0.17± 
0.02. This profile shows this self-similar behavior to the 
resolution limit of the simulations at 100 kpc/h. 

We choose conservative models the of dark matter halo 
profile in calculating Jaq. that are consistent with CDM 
in addition to observational constraints of the MW halo. 
The MW CDM halo profile has been constrained with dy- 
namical measurements by Klypin et al. [H] and Battaglia 
et al. [H] using NFW-type halo profiles. Recent results 
by Catena & Ullio (CU) [10] have applied such dynami- 
cal constraints to NFW as well as Einasto profiles of the 
MW with a Markov-Chain Monte Carlo parameter ex- 
ploration of the degeneracies in fit parameters in these 
models, providing the local (solar) density estimates of 
the MW halo as well as the global profiles for these halos. 
The results of CU allow using not only minimal models 
in MW halo mass, but variations in other parameters as 
well. 

There exists an overwhelming amount of results from 
verified numerical simulations that halos in CDM cos- 
mologies generally form cusped profiles of the form ex- 
emplified by the NFW and Einasto profiles [TTJ H7P0]. 
The original ansatz of CDM has been shown to be consis- 
tent with a large range of scales of observations, includ- 
ing the cosmic microwave background [50] . large scale 
clustering of galaxies and dark matter, ^100 Mpc/ft,, in 
the Sloan Digital Sky Survey [5T], to the small scales of 
the intergalactic medium measured in the Lyman-a for- 
est, ~50 kpc/h [52H54] , We do not consider isothermal 
and Burkert profiles, since cored profiles of these forms 
never arise in CDM halo formation simulations of halos 
of the mass of the MW. We also do not consider the 
more shallow-cusped "Kravtsov" profiles originally from 
Ref. |55j since following work by the same authors found 
sufficient force resolution but insufficient mass resolution 
in the original simulations, leading to the spurious shal- 
lower profiles |56| . 



Galactic Center Observations 



The lower and upper limit of the integration cc m i n , x max 
are set by whether it's the MW or an extragalactic halo. 
The limits chosen are somewhat arbitrary since galactic 
halos are embedded and continuous with a network of 
lower-density filamentary structures. However, the con- 
tribution beyond the scale radius is minimal, and we sim- 
ply choose a cutoff at the nominal virial radius for the 
MW. The average of the integrated mass density is 



J An 



(6) 



over the relevant sky region, where Af2 is the angular size 
of the observation region. Although the NFW profile is 
divergent at the GC, the average over a finite observation 
region is insensitive to the exact nature of the inner cusp 
for the 7-ray observations here. We regularize the NFW 
profile by fixing a constant density of the profile within 
''halo < 10~ 10 kpc to the density at p(rh a io = 10~ 10 kpc) 
in Eq. [6j which has negligible effect on the integrated 
average. 

One of the most stringent limits on dark matter anni- 
hilation arise from observations of the MW Galactic cen- 
ter ridge (GR) by the High Energy Stereoscopic System 
(HESS) telescope [57], an array of atmospheric Cerenkov 
telescopes in operation in Namibia 6 [Ml HU EH [32j l58l - 
150] . The HESS GR observations are for an angular re- 
gion -0.8° < t < 0.8° and -0.3° < b < 0.3°. We inte- 
grate Eq. ([6| directly, through a Monte-Carlo integration 
method over the nontrivial geometry. Note that in HESS 
GR observation subtracted a background from a region 
-0.8° < I < 0.8° and 0.8° < b < 1.5° to remove cos- 
mic ray backgrounds. This also removes any concurrent 
dark matter signal in the background region, which we 
incorporate appropriately by integrating the signal from 
the background subtracted region —0.8° < t < 0.8° and 
0.8° < b < 1.5° numerically, also through a Monte Carlo 
integration method for the arbitrary geometry. 7 The full 
Monte Carlo integration particularly affects the Einasto 
profile case due to the shallower Einasto profile which 
falls off more slowly than NFW in the background region, 
and therefore increases the background subtraction. 

We calculate extremal cases for Jc^i hgr towards the 



For comparison, we include constraints from observa- 
tions towards the center of our Galaxy. The integrated 
mass density squared along line-of-sight x towards Galac- 
tic coordinates (6, t) from the GC is 



j{b,e) = h 



p 2 {r m \{b,i,x))dx, 



(4) 



where J = 1/ 8.5 kpc (0.3 GeV cm 3 )~ is a normal- 
ization that makes J unitless and cancels in final expres- 
sions for observables, and 



r ga i(b 7 £,x) — J Rq — 2xRq cos(£) cos(6) 



(5) 



6 http://www.mpi-hd.mpg.de/hfm/HESS/HESS.html 

7 Note that we do not adopt an often-employed circular area 
approximation for the HESS GR observational geometry, or a 
point approximation for the background subtraction, e.g. in 
1 131 1141 . The circular approximation underestimates the signal 
by approximately a factor of two, primarily because the obser- 
vational area is overestimated while the signal is centrally con- 
centrated. The circular area approximation with angle ?/> = 0.8° 
has Afi circ = 2tt(1 - cosi/i) = 6.1 X 10~ 4 sr while the HESS 
GR observation was over an area ASljjESS = 2.9 X 10 -4 sr. The 
background approximation of a point subtraction at ip = 0.8° 
overestimates the background integral, yet this is partially offset 
by the smaller area of the background Af2pjESS.bg = 3.4x 10 -4 sr 
than in the point approximation. 
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HESS GR field of view, such that 

555 CU Einasto Minimum 

1050 Battaglia NFW Low Mass 

Jao,hgr = { 2500 CU NFW Low c (7) 

13100 Battaglia NFW High Mass 

„ 14700 CU NFW High c. 

For the minimal Einasto case, we use a softer conserva- 
tive value of the log-slope that is consistent with the fits 
of CU, with a E = 0.22, r s = 21 kpc, r© = 8.28 kpc and 
Pq = 0.385 GeV cm -3 which are the values we adopt 
as our conservative minimal case of the MW halo profile 
in the HESS GR field of view. Our adopted maximally 
stringent case is the NFW profile from CU which has the 
highest concentration, c. High concentration enhances 
the signal toward the GC, yet is anti-correlated with halo 
mass therefore corresponding to a low mass MW 
halo, with c = 24.6 and M vir = 1.23 x 10 12 M . The CU 
NFW low c case has c = 13.9 and M vir = 1.86 x 10 12 M Q 
(we do not employ this model in our constraints). Both 
CU NFW profiles adopt r© = 8.28 kpc. The high and 
low mass NFW cases from Battaglia et al. [45] corre- 
spond to the 68% CL limits of the MW halo model: 
M vlr = 0.6 x 10 12 Mq (low mass), M vir = 2.0 x 10 12 M Q 
(high mass). Both of these cases have other parame- 
ters fixed at i? vir = 255kpc, c = R v - lr /r s — 18, and 
r Q = 8 kpc. The Battaglia et al. NFW extremal profiles 
and corresponding GC J are consistent with those from 
CU, see Eq. Q. The CU and Battaglia et al. NFW pro- 
files are also consistent with and comparable to the high 
and low mass MW halo models of Klypin et al. [33] . In or- 
der to allow for comparisons with previous work, we also 
evaluate the constraints for the HESS GR for the case 
of the canonical NFW halo with p Q = 0.3 GeV cm~ 3 , 
i? = 8.5 kpc, r s — 20 kpc, where J7aq,hgr = 1620, and 
the constraint in parameter space is plotted as a dotted 
line in Fig. [2] 

Our results are less stringent for the Einasto case com- 
pared to Refs. [231 [35] because the signal used in that 
work did not remove the signal from the background sub- 
traction region. Note also that our minimal case of a "CU 
Einasto Minimum" is significantly more stringent in its 
constraints than the adopted values in, e.g. Refs. [T3l[Ti] 
which adopted inaccurate MW profiles from Ref. [55] 
which had known insufficiencies in numerical mass res- 
olution as discussed above. We do not include here con- 
straints from the smaller region HESS GC observations 
presented in Ref. |61j since they are weaker than the 
HESS GR observations for nearly all channels and re- 
gions in parameter space |32j . 



D. The Dark Matter Signal in the Isotropic Diffuse 
Background 

The MW halo itself contributes to the isotropic back- 
ground in an irreducible fashion due to the smooth and 



unresolved portion of the MW halo. We use a very con- 
servative limit in the required contribution of the MW 
Galactic halo to the isotropic diffuse background. Specif- 
ically, we assume that the entire diffuse MW signal is 
at minimum the value opposite the GC from a smooth 
halo, with a minimal substructure boost. We take that 
Jin is given by the anti-GC ^7(0,180°) averaged over 
the full diffuse observation region Af2, J\ so ~ 0.62. The 
signal towards the most Galactic-centric portion of the 
Fermi-LAT diffuse observation region is over two or- 
ders of magnitude larger than that away from the GC, 
4/(10°, 0) ~ 78. Therefore, our approach is extremely 
conservative, and leads to a robust bound since the re- 
moval of a spatial dependence in the Fermi-LAT diffuse 
spectrum may remove spatially-dependent components 
of the predicted signal. 

A cold dark matter halo necessarily has enhanced an- 
nihilation from substructure in the dark matter halo. 
Early work had found MW substructure to contribute to 
boost factors of the total MW annihilation signal at over 
two orders of magnitude greater than the smooth halo 
contribution |62j . More recent work from the Aquarius 
simulations find similar boost factors at 232 times the 
smooth contribution [63] . The minimal level of enhance- 
ment from numerically resolved substructure has been 
found to be a factor of ^2 from the Via Lactea and Via 
Lactea II simulations [531 [55] ■ However, substructure is 
expected to exist in CDM halos at potentially over twenty 
orders of magnitude in mass scale below a MW parent 
halo. Therefore, we go beyond the approximations from 
resolution limitations. 

In order to provide a conservative estimate of the en- 
hancement of annihilation due to unresolved substructure 
in the Fermi-LAT observation, we employ the analytic 
density probability distribution function (PDF) proposed 
by Kamionkowski & Koushiappas [66 . This method was 
recently applied to the enhancement of annihilation sig- 
nals in MW-type halos with a calibration of the PDF 
with the Via Lactea II (VL-II) simulations [57]. That 
work found typical boost factors for the full MW halo 
are approximately 17. 

In order to quantify the lower limit of the boost factor, 
we extremize the PDF to the lowest values consistent 
with the CDM halo calibration. The mean halo density 
from the PDF is given by 



p(r) 



P P(p)dp 



(8) 



10 

fsPh + 

(1 - fs)Ph 



1+q 



a = 0, 



(l-/> h ln(W); a^0. 



(9) 



The fraction f s of the halo volume is filled with a smooth 
dark matter component with density ph- The power law 
tail is calibrated to simulations to be 



(1 - / s ) = 7 x 10" 



p{r) 



p{r = 100 kpc) 



-0.26 



(10) 
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FIG. 3: Shown is a comparison of our work to the Fermi-LAT 
collaboration, Abdo et al. [12] constraints from the extra- 
galactic signal in the diffuse extragalactic background from 
dark matter annihilation into the bb final state. Our limits 
from the extragalactic background with not boost enhance- 
ment of annihilation are comparable to and consistent with 
that of the Abdo et al. "MS-II Res" model due to similar cho- 
sen halo mass cut-off values. Our total 95% CL limit shown 
here and used in Fig.[2]comes from Galactic and extragalactic 
contributions, with a conservative boost factor, as described 
in the text. Shown also is our 95% extragalactic limit alone, 
and the conservative model in Abdo et al. of "MS-II subl." 



for radii greater than 20 kpc, where the overwhelmingly 
dominant portion of the substructure boost arises. The 
normalization of this tail is lower at smaller radii, < 
20 kpc, but the contribution to the cumulative boost from 
these radii is minimal. Simulations find that a — 0±0.1. 

The definition of the annihilation boost factor includ- 
ing the CDM density PDF, as a function of galactic ra- 
dius is 



B(r) = 



Jp 2 dV 
J[p(r f dV ' Jo 



P{p,r) 



rdp. (11) 



Explicitly, for this PDF, the boost factor is 
B(r) = f s e A2 



Hl-fs) 



1 



1 - a 



Ph 



(12) 



A 2 

The first term f s e is due to the variation in the smooth 
component, and since A < 0.2 from simulations, it con- 
tributes to the overall boost factor by only a few percent 
and we ignore it. The second term is the boost factor 
due to substructure. The total boost in the galactic halo 



B(< R) = 



J B(r)p(r) 2 r 2 dr 
Jo* p{r) 2 r 2 dr 



(13) 



We take the extremal case to be where the slope of the 
power-law tail is steeper than the mean found in simu- 
lations, a — 0.1. We take p max to be p s of the earliest 
forming halos, which have low concentrations c ~ 2 — 5, 
which for centrality of this concentration distribution we 
take c = 3.5, forming at z c — 40. Note that using p s is 
extremely conservative since the smallest subhalos likely 
have cuspy profiles themselves that reach densities much 



higher than their scale density p s . The integral, Eq. (13) 



must be evaluated numerically, and for our adopted con- 
servative case is 



B{< i? vir ) = 6.6, 



(14) 



for i? v i r = 255 kpc. Note that this total boost is signif- 
icantly lower than typically found in analytic extrapola- 
tions of MW halos which find boost factors of over two or- 
ders of magnitude. We employ this minimal boost factor, 
Eq. (14 1, in our signal calculations. It should be noted 



that the boost in the local (solar) region is minimal, on 
average, but including the boost due to substructure can 
affect and improve the modeling of leptonic cosmic-ray 
annihilation signals (67] [68] . 

A similar modeling to the one we adopt here of 
the density PDF calibrated to cosmological volume 
(100 Mpc/h) 3 simulations found larger scatter of the 
power law tail at mass scales approaching MW scale halos 
with boost factors between 2 and 2000 [55] • The mass res- 
olution of the Millenium Simulation II (MS-II) employed 
in Ref. [69] (particle mass 6.89 x 10 6 Mq/K) is approxi- 
mately 2400 times less than that in VL-II (particle mass 
4100 M ), and the MS-II results did not extend to PDF 
resolution of MW size halos. This may contribute to the 
increased PDF scatter for lower mass halos, and there- 
fore we adopt the VL-II calibrated PDF of Ref. [57]. Note 
that a complete systematic analysis of the lower limit of 
the substructure PDF contribution to the boost has not 
been performed since some of the earliest work on sub- 
structure boost factors, which found total boost scatter 
at the level of 40% 1 62 . Such a lower scatter is consistent 



with what we find, Eq. (14) 



The cosmological contribution to the isotropic diffuse 
7-ray background from annihilating dark matter has been 
studied for some time [701 [7T] . The diffuse flux from 
annihilation in cosmological extragalactic halos is 



~d~E 



WAV) 



(/DM^n 



AS! b s -ffo 



> f{z){l + zYdN{E>) UE , 

h(z) dE> 6 ' ( 5) 



where Hq = 70 km s~ x Mpc~ x is the Hubble constant, 
Q m is the matter density in units of the critical den- 
sity, p cr it, and the fraction of matter in dark matter 
is /dm = ^dm/^dm + fife) ~ 0.833, where the frac- 
tion of critical density of the dark matter we take is 
Hdm = 0.237, and baryon density f2& = 0.0456 [50] . 
The uncertainties in the halo modeling are much larger 
than the errors on the cosmological parameters we adopt 
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FIG. 4: Shown are constraints on regions of parameter space in annihilating dark matter models consistent with interpretations 
of the PAMELA positron excess (light pink region in all panels), and feature in the Fermi-LAT e + /e~ spectrum (red region 
in fi + fi~ and t + t~ panels), from Ref. [24] (all 99% CL). All MW constraints and signals are for Einasto profiles, and at 95% 
CL. Our analysis of the Fermi-LAT isotropic diffuse background is shown along with our reanalysis of constraints from HESS 
observations of the Galactic Ridge (labeled HESS GR). We show several other recent constraints for comparison: ICS radiation 
constrained by Fermi-LAT data from the 3° x 3° Galactic center exclude the region within the triple-dot-dashed regions in all 
panels ,|T; Galactic radio synchrotron observations exclude the region within the dotted line, in all panels in the central 
fi + fj,~ panel, we show the exclusion regions final state radiation constraints in Fermi-LAT observations of Draco (dashed line), 
and ICS radiation constraints from Fermi-LAT observations of Ursa Minor (dot-dashed line) [7|. All PAMELA models with 
m x < 1 TeV are firmly excluded by the lack of a > 20% drop in the Fermi e + /e~ spectrum below 1 TeV |24| . 



here. However, much previous work sets /dm = 1, which 
overestimates the extragalactic annihilation signal by ap- 
proximately 30%. We take a flat universe, with 
and h(z) = [(1 + zfn m + f2 A ] 1/2 . The factor e - r ( 2 < B '> 
takes into account attenuation of 7-rays along cosmolog- 
ical distances, for which we use the results of Gilmore 
et al. [72]. The factor f(z) accounts for the increase in 
density squared during halo growth and the redshift evo- 
lution of the halo mass function. We adopt a fit to this 
evolution such that [7T1 1751 174] 



f( z ) = y io - 9 [ ox P(- - 9z )- 1 ]-°- lfo . 



(16) 



The halo internal density profile sets /o, and for the 
Einasto case, we find fo — 3 x 10 4 . We also include a the 



boost factor here of Eq. ( 14 ) to the full extragalactic com- 



ponent, consistent with our adoption for the MW. The 
extragalactic contribution to the diffuse isotropic back- 
ground is subdominant to the Galactic contribution, at 
the level of 10%-20% of the total signal for a given model. 
For comparison, shown in Fig. [3j our model of the diffuse 
extragalactic background is similar to that named "MS- 
II Res" in Ref. [TJ], due to a comparable low mass halo 
cutoff, excluding our adopted boost factor. 

We redshift the energy distribution of the annihilation 
photons as E = E' / (1 +z), where E' is the energy of the 
photons at the cosmological source. Eq. ( [To] ) uses a min- 
imal limit of the halo mass of 10 6 M /h, a conservative 
lower limit that only uses cosmological halos known to fit 
cosmological determinations of the halo mass functions 
with no extrapolation [75]. Note that the limit from the 
diffuse flux greatly increases with an extrapolation to ex- 
tremely low masses such as 10 -9 M Q employed by Hiitsi 



ct al. 



With our framework, we use the observed cos- 



mological diffuse spectrum as the maximum amount that 
the annihilation signal can be an a given energy with 95% 
CL limits to the full spectrum for any given m x . For the 
95% CL limit, we sum take the bin-summed x 2 < 1-282 
for this form of a one-sided upper limit. These limits are 
shown in Fig. [5] 



E. PAMELA and Fermi 
Electron/Positron-Spectrum Motivated Models 

The PAMELA observation of an increase on the 
positron fraction at 10 to 100 GeV, in combination with 
a feature in the shape of the Fermi e + /e _ spectrum at 
~1 TeV could be consistent with the production of the 
high energy e + /e~ in the products of dark matter anni- 
hilation. Since annihilation modes to charged particles 
also produce photons from bremsstrahlung, the signal is 
should also be observed or constrained by the Fermi-LAT 
isotropic diffuse spectrum. Channels through the quarks 
and massive W and Z bosons are strongly excluded by 
a number of observations |24j . We examine here con- 
straints on two-body charged annihilation modes from 
the Fermi-LAT isotropic diffuse background, with 95% 
CL exclusions shown in Fig. |2J The PAMELA 99% CL 
preferred region is in pink, and combined Fermi-e + /e _ 
preferred 99% CL region in red, from Ref. |24) . 

We again choose the most conservative cases for the 
MW halo profile, as in §11 D] with a shallow Einasto pro- 
file as that incorporated in Eq. (JtJ) , and high mass cut- 
off of the halo mass function at M m - m = 10 6 Mq. The 
Fermi-LAT diffuse spectrum excludes a large part of the 
parameter space for the PAMELA and Fermi-e + /e~ 99% 
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m x [GeV] 

FIG. 5: Interpretations of PAMELA and Fermi e + /e~ with an 
intermediate dark force carrying particle cj> allowing for dark 
matter annihilation into four lepton final states. The upper 
(lower) two panels are 4e (4fi) final states with — 0.3 GeV 
and rrirf, = 0.8GeV. Galactic radio synchrotron observations 
exclude the region within the dotted line, in all panels [24] , 
All PAMELA models with m x < 1 TeV are firmly excluded 
by the lack of a > 20% drop in the Fermi e + /e~ spectrum 
below 1 TeV [Hj. 



CL regions, and when combined with the HESS Galac- 
tic Ridge constraints, exclude all interpretations of the 
PAMELA positron fraction and Fermi-e + /e~ feature as 
arising from dark matter annihilation into two-body stan- 
dard model particle final states. 

We also consider four-lepton annihilation modes that 
could occur through an intermediate force carrying vector 
or scalar boson, tj>, with mass < 1 GeV. We take 
scalar bosons of two cases, = 0.3 GeV and = 
0.8GeV, into either 4e or 4/z final state modes. They are 
shown in Fig. [5j For the more massive <j> particle, — 
0.8GeV, the resultant lepton states more energetic and 
therefore more constrained by their final state radiation 
by the Fermi-LAT diffuse and HESS GR observations. 
For either </> mass, there exist regions of the parameter 
space that are not excluded at the 95% CL level by the 
gamma-ray data we consider here. 

For reference, we also show several other constraints 
that have been placed in the literature on these models 
also using Einasto-type profiles for the MW halo. For 
the two lepton modes, we show the region excluded by 
ICS radiation of the original final state radiation prod- 
ucts scattering off of cosmic microwave background pho- 
tons, as constrained by Fermi-LAT data from the Galac- 
tic poles (triple-dot dashed) , from Ref. [9] , and originally 
described in Ref. [75]. We also show the region excluded 
by synchrotron emission not observed in Galactic radio 
observations (dotted line) from Ref. [24 , originally de- 
scribed in Ref. [32] ; for the case of x + X ^ + + M~ > 
we show the exclusion regions from final state radiation 
constraints in Fermi-LAT observations of Draco (dashed 
line), and ICS radiation constraints from Fermi-LAT ob- 
servations of Ursa Minor (dot-dashed line) [Jj. In the 
four lepton modes, we show Galactic radio synchrotron 
constraints [23]. Also, it should be noted that all mod- 
els with m x < 1 TeV are firmly excluded by the lack 
of a > 20% drop in the Fermi e + /e~ spectrum below 
1 TeV [23] • The ICS constraints from public data of 
Fermi-LAT observations towards the Galactic poles [9 
are comparable to and, in many regions, stronger than 
our constraints, though such constraints do depend on 
ICS modeling of the signal in these regions. 



III. CONCLUSIONS 

The LAT aboard the Fermi Gamma-Ray Space Tele- 
scope is opening a new window to detecting the nature 
of dark matter or constraining it. We have derived con- 
servative constraints on WIMP dark matter annihila- 
tion cross-sections for all standard model channels aris- 
ing from the observed isotropic diffuse background by the 
Fermi-LAT. We have also re-examined constraints from 
the HESS Galactic ridge observation. In contrast with 
previous work, our methods use very conservative mod- 
els for the cosmological and halo structure giving rise to 
the signal, while still remaining consistent with structure 
formation in a cold dark matter framework of a WIMP 
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candidate. We use such conservative models in the dark 
matter halo profile in the Galactic signal contribution, 
and conservative lower-limits to the surviving halo mass 
function in the Galactic and extragalactic contributions 
to the isotropic diffuse background. For example, we do 
not consider isothermal halo profiles which do not arise 
in any CDM halo formation calculations. 

Our conservative models are therefore stringent and 
robust constraints on any WIMP-like annihilating dark 
matter particle candidates. The constraints are at or near 
the thermal cross-section for annihilation to hadronic 
modes, quarks and gluons. For low masses, 5 GeV < 
m x ~ 10 GeV, the constraints exclude the expected 
thermal-annihilation scale (ctav) ~ 3 x 10~ 26 cm 3 s _1 
into hadronic final states. Such low masses are not a 
priori excluded by previous particle physics or cosmolog- 
ical constraints [771 [78]. Therefore, such constraints are 
approaching an observational lower bound on the dark 
matter mass for the thermal-relic WIMP. The constraints 
from the isotropic diffuse background presented here are 
more stringent than the constraints from the diffuse back- 
ground limits of the Fermi-LAT collaboration when in- 
cluding only the extragalactic component [12] . and com- 
parable to Fermi observations towards dwarf galaxies [7J , 



galaxy clusters [79], and the Galactic poles [5], while be- 
ing more stringent than full-sky photon limit methods [8] . 
With increased integration time over the life of the Fermi 
Gamma-Ray Space Telescope and further constraints on 
dark matter structure in local and extragalactic sources, 
the isotropic diffuse spectrum or other Fermi-LAT obser- 
vations will either reveal the 7-ray products of WIMP- 
like dark matter annihilation, or exclude this class of can- 
didates. 
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